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Abstract. Polar stratospheric clouds (PSCs) and cirrus have been investigated using thermal 
emission measurements at 10.8 and 12 urn wavelength (channels 4 and 5) from the Advanced 
Very High Resolution Radiometer (AVHRR). The AVHRR signal was evaluated from a 
theoretical basis to understand the emission from clear and cloudy skies, and models were 
developed to simulate the AVHRR signal. Signal simulations revealed that nitric acid PSCs are 
invisible to AVHRR, while ice PSCs and cirrus are readily detectable. Methods were developed 
to retrieve cloud optical depths, average temperatures, average effective radii, and ice water 
paths, from AVHRR channels 4 and 5. Properties of ice PSCs retrieved from AVHRR were 
compared to values derived from coincident radiosondes and from the Polar Ozone and Aerosol 
Measurement II instrument, showing good agreement. 


1. Introduction 


^ This work explores the detection of polar stratospheric clouds (PSCs) and cirrus using 
' measurements from the Advanced Very High Resolution Radiometer (AVHRR). Interests in cirrus 
include their role in radiative forcing of the climate system and their role in the chemical balance 
of the atmosphere [e.g., Sassert et al., 1989; Borrmann et al., 1996]. Because cirrus usually occur 
below the spontaneous freezing temperature of water, they are typically composed of ice crystals 
Jensen etal., 1994]. PSCs have been implicated as key factors in the loss of ozone under cold polar 
conditions. Observational and theoretical evidence suggests that PSCs can be composed of either 
liquid ternary H 2 SO 4 -H 2 O-HNO 3 aerosols (LTA), nitric acid trihydrate (NAT), or ice [e.g., 
Tolbert, 1996; Carslaw et al., 1994]. 
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AVHRR measures reflected and/or emitted radiation at nadir (±50°) in five spectral bands 
between 0.58 and 12.5 fun wavelength (X). A series of AVHRR sensors on board polar-orbiting 
satellites have been collecting data since 1979, and provide excellent temporal and spatial coverage 
of the polar regions. Poleward of 60°, each point on the earth’s surface is sampled in two or more 
consecutive passes twice a day for each of the two polar-orbiting satellites. These images can be 
combined to provide daily coverage of each polar region with horizontal resolution of 1.1 km at 
nadir. The spatial overlap of consecutive images allows for some tracking of cloud movement, 
although for a reladvely short time period when compared to geostationary satellite data. AVHRR 
measurement characteristics, including the ability to measure with or without sunlight, offer unique 
advantages for PSC detection. Garcia, et al [1995] and Pagan [1996] have shown that AVHRR 
measurements of 12 pm brightness temperature can be used to detect optically thick PSCs when 
underlying Earth surface temperatures are well above PSC temperatures. Most PSCs, however, 
are optically thin and occur over a variety of underlying surface conditions, and cannot be detected 
using a single AVHRR channel. 

Cures and ice PSCs can have similar characteristics, so it is reasonable to expect that 
techniques established for cirrus characterization could be applied to PSC measurements. Inoue 
[1985, 1987] found that cures clouds can be identified over tropical oceans using the difference 
between channel 4 and 5 brightness temperatures (T 4 and T 5 ). This brightness temperature 
difference (BTD = T 4 - T 5 ) was then used by Yamanouchi et al. [1987] to detect thin cirrus over 
Antarctica. The radiative contrast between channels 4 and 5 has been investigated as a tool for 
inferring physical cloud properties such as temperature and particle effective radius [e.g., Wu, 1987; 
Parol et al., 1991; Lubin and Harper, 1996]. Parol et al. concluded that while cloud effective radii 
cannot be accurately determined, that useful limits can be placed on its range. Lubin and Harper 
inferred the temperature, ice water path, and effective radius of Antarctic cirrus, and suggest that 
these properties can be determined with some certainty for conditions relevant to cirrus and PSCs. 

For this work, the AVHRR channel 4 and 5 radiances were evaluated from a theoretical basis 
to understand the emission from clear and cloudy skies. Methods were developed for retrieving 
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cloud extinction optical depths, average temperature (T c ), average effective radius (r c ), and ice 
water path (TWP), from the channel 4 and 5 radiances. These retrievals were applied to AVHRR 
cloud measurements for comparison to values derived from radiosonde measurements and to 
measurements from the Polar Ozone and Aerosol Measurement (POAM) II instrument. 

2. Radiative Transfer and Signal Simulation 

The nadir emission from a column of atmosphere containing cloud, L m (X), can be expressed 
in terms of net radiances (L) and transmissivities (t) 

Ljn(^) = { [L E (X) T t (X) + L At (X)] x c (k) + Lc(X) } t u (X) + L Au (X) (1) 

where L E (X) is due to the Earth’s surface, T t (X) and L At (X) are due to clear air below the cloud layer, 
x c (X) and Lc(X) are due to the cloud layer, and t u (X) and L Au (X) are due to clear air above the cloud. 
Each L and x is an integral in height and wavelength, and the expanded relationship is tedious 
compared to (1). The terms in (1) are now explored. 

\ 

2.1. Earth’s Surface 

Radiance emitted from the Earth’s surface can be expressed in terms of the surface emissivity 
E 0 (X) and the blackbody radiance (L*) at surface temperature T E , L E (X) = L*(T E , X) eo(X). For 
typical surface cover underlying PSCs (sea water, snow, and ice), values of Eo(X) are close to unity 
and are roughly constant in the 10 to 12 pm region (Table 1). Thus, the spectral content introduced 
to L m (X) by L e (X) is primarily that of the Planck function. 

X 

2.2 Clear Skies 

Atmospheric absorption and emission were determined using rigorous line-by— line 
calculations [Gordley et al., 1994] for interfering gases in the channel 4 and 5 bandpasses. In order 
of radiative importance, the species considered were H 2 0, C0 2 , 0 3 , HN0 3 , CC1 3 F, and CC1 2 F 2 . 
Frequency integrals were convolved over the AVHRR spectral response functions to more 
accurately simulate an AVHRR measurement. The vertical integrals considered transfer from the 
Earth’s surface to altitudes where gaseous contributions become negligible (about 50 km). Clear 
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sky radiances were simulated for conditions ranging from polar winter to tropical. For this purpose, 
vertical profiles of temperature, pressure, and gas mixing ratios were taken from climatologies. 
These calculations show that the atmospheric contribution is determined mainly by H 2 O and CO 2 
in the lower troposphere, and that transfer above the tropopause can be safely neglected (i.e., x u (X) 

“ 1 311(1 L AuM “ °)- Thus, the clear sky emission, L CS (X), is determined by the Earth’s surface 
and radiative transfer in the troposphere, L CS (X) = L E (X) x t (X) + L At (X). 

The clear sky results allow some simplifications of (1) which lead to a more general 
expression for cloudy sky emission. L CS (X) could be measured from clear sky adjacent to a cloud, 
allowing the cloudy sky radiance to be expressed in terms of two measurements and two cloud terms. 
In a more general sense, the radiance illuminating a cloud from beneath may be affected by 
uncharacterized cloud layers at lower altitude. Considering that this background radiance, L bk (X), 
can be characterized by AVHRR measurements away from a cloud of interest, the cloudy sky 
emission (1) is reduced to 

LmM = L bk (X) x c (X) + Lc(X) (2) 

where L bk (X) can be from a clear background, L bk (X)=L cs (X), or cloudy background. The form of 
(2) is similar to expressions derived by Parol et al. [1991] and Lubin and Harper [1996] for cloudy 
sky emission. Figure 1 depicts the concept of clear and cloudy sky signals adopted for this work. 
Potential limitations to this approach come quickly to mind. However, (2) eliminates the need to 

accurately characterize the Earth’s surface and atmosphere for input to rigorous radiative transfer 
calculations. 

The dependence of clear sky radiance on telescope scan angle should be understood since 
values of L CS (X) and L^X) could be taken from measurements at different scan angles. Clear sky 
calculations for a range of scan angles considering tropical and polar conditions show that sensitivity 
to scan angle increases for optically thick (e.g., tropical) atmospheres. Transmissivity is noticeably 
reduced at large scan angles. However, emission increases in proportion to x t (X) so that net radiance 
is relatively insensitive to scan angle (Figure 2). 
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Figure 1. A conceptual model of thermal emission from clear sky and one or two cloud layers 
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Table 1. 

Various Properties Relevant to AVHRR Channels 4 and 5. 

Channel 

Bandpass Wavelengths (pm): 
Lower Limit, Center, Upper Limit 

Cloud Particle: Complex 
Refractive Index 

Earth Surface: 
Emissivity 

4 

10.3, 10.8, 11.3 r 

ice: 1.09 +0.1 8i (W) 
NAT: 1.63 +0.13i (T) 
LTA: 1.58 + 0.24i (H)* 

snow: 0.993 (M) 
sea ice: 0.985 (G) 
sea water: 0.99 1 (G) 

5 

11.5, 12.0, 12.5 r 

ice: 1.28+ 0.4 li (W) 
NAT: 1.53 + 0.27i (T) 
LTA: 1.59 + 0.27i (H)* 

snow: 0.972 (M) 
sea ice: 0.970 (G) 
sea water: 0.986 (G) 

The information was taken from Gao et al. [1994] (G), MOSDIS UCSB 
www.icess.ucsb.edu (M), Warren [1984] (W), Hervig etal. [1997] (H 
[1994] (T). Refractive indices and emissivities are at the band center v 
*LTA at 192 K composed of 10 wl% H 2 S0 4 - 35 wt.% HN0 3 . 

'•'Values for the NOAA-12 sensor. 

Emissivity Library, 
), and Toon et al. 
/avelengths. 


23* Cloud Optical Properties 


The extinction coefficient (P e ) for a population of particles is calculated by integrating the 
extinction efficiency (Q e s optical cross section / geometric cross section) over the size distribution: 


m = 


I 


Kr2 Qe(r,X,m) n(r ) dr 


(3) 


where m is particle refractive index (Table 1), n(r) is the particle size distribution, and rj and r 2 
are the smallest and largest radii in the size distribution. The scattering coefficient (f3 s ) is calculated 

by substituting the scattering efficiency (Qs) into (3), and the absorption coefficient (0J is the 
difference (3 e ~p s . 


The cloud particle size distribution can be described using a lognormal distribution for PSCs 
[Deshler et al., 1994] or for cirrus [Jensen et al., 1994]. While real cirrus or PSC size distributions 
can be quite different from lognormal, the lognormal form is useful for considering a broad range 
of conditions in modeling cloud radiative properties. The lognormal distribution is described by 


nlrfdr -f N ... 

where n(r) is the differential concentration in the radius interval dr, k is the number of modes in the 
distribution, N; is the number concentration, wj is the distribution width, and r, is the median radius. 
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for each of i modes. Cloud volume density (V, particle volume per unit volume of air) and effective 
radius (r e , surface area weighted mean radius) can be derived from the size distribution. These 
important quantities are introduced here for convenience. 

r e = r exp[5 (ln 2 w)/2 ] ( 5 ) 

V = N (4/3) x P exp(8 (In 2 w)/2 ] (6) 

In this work, optical cross sections for (3) were determined using Mie theory, which explains 
the scattering and absorption by spheres. While cirrus and PSC particles are probably nonspherical, 
Mie theory can still provide reasonable estimates [e.g., Ebert and Curry, 1992]. To test this 
assumption, calculations based on Mie theory were compared to results from Fu and Liou [1993] 
(FL) who used ray tracing and spheroidal solutions to calculate extinction cross sections for 
randomly oriented hexagonal ice crystals. FL computed integral extinctions, 0 C (X), for 1 1 measured 
size distributions and fit the ratio f3 e (X)/V versus effective radius for numerous broad spectral bands 
across the visible and infrared. The bands used by FL did not match AVHRR channels 4 or 5, so 
one of their bands between 10.2 and 12.5 pm (central X = 11.3 pm) was chosen for the comparisons. 
Using the Mie-Iognormal model, P e (11.3) and V were calculated for a range of lognormal 
distributions. The resulting volumes encompass the FL results (Figure 3a), so it should be possible 
to mimic their size distributions by carefully selecting lognormal parameters. Because volume and 
extinction are roughly proportional, we expect the ratio p e (X)/Vto be somewhat independent of size 
distribution. To the extent that this is true, the comparison of 0 e (l 1.3)/V is a good measure of the 
differences between Mie extinction and the more rigorous calculations of FL (Figure 3b). The Mie 
extinctions were scaled by the ratio of volume from FL over the lognormal volume, to account for 
size distribution differences. The scaled Mie extinctions are 10 to 30% greater than FL (Figure 3c). 
However, it is evident in Figure 3c that the Mie-lognormal model could exactly reproduce FL by 
adjusting the choice of lognormal size distributions. Since FL considered only 11 size distributions, 
there is no reason to expect that their results are more versatile than results based on a large family 

of lognormal size distributions. Thus, this work proceeds using the readily available 
Mie-lognormal model. 
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10 100 
effective radius (pm ) 


Figure 3. Results from Fu and Liou 
[1993] (FL) compared to calculations 
based on Mie theory and lognormal size 
distributions (ML). Both FL and ML 
used ice refractive indices from Warren 
[1984]. In Figures 3a-3b, the groups of 
dots represent results from a range of 
lognormal distributions (w = 1.5, 2, and 
2.5; r = 0.1 to 87.1 pm), and dashed 
lines are fits to these results. The 
differences are between FL and the ML 
fits (dashed curves), a) Volume density 
versus effective radius from the 
lognormal distributions compared to 
FL. b) Comparison of extinction over 
volume. c) Comparison of the 
extinction coefficient. The fitted ML 
curve shown here was scaled by the rado 
of FL volume over lognormal volume 
(solid and dashed curves in Figure 3 a). 


Cloud transmissivity is determined by radiadve loss due to scattering and absorprion 


(extinction) in the cloud layer. 


x c (X) = exp(-(3 c (X) Az) (7) 

where Az is geometric cloud thickness. Cloud emissivity is determined by a cloud’s absorptive 
properties. 


e c (X) = 1 - exp(-f} a (X) Az) (8) 

and cloud radiance is determined by the cloud temperature (T c ) and emissivity, > 

Lc(K) = L*(T C) X) e c (X) (9) 

where L*(T C , X) is the blackbody radiance. Optical depth, o e (X) = 0 c (X)Az or a a (X) = 0 a (X)Az, can 
be substituted into the above relationships. The cloudy sky emission (2) can be restated in terms of 
a single optical depth by defining the ratio /^ e (X) = a a (X)/o e (X) and substituting (7), (8), and (9). 

Lm (X) = 4*(X) exp[-o e (X)] + L*(T C , X) { 1 - exp[-o c (X) (X)] } (10) 

If Rac(X) can be treated as constant, (10) defines a closed system, relating two measurements and 
two unknowns. Calculated values of RaefX) are bounded between 0.45 and 1.0 for ice refractive 
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indices and a range of size distributions representative of PSCs or cirrus (Figure 4a). It can also be 
shown that R ae ( 10 . 8 ) and R ae (12) are strongly correlated for ice indices and a wide range of size 
distributions (Figure 4b). While fixing R ac (X) may or may not prove acceptable, its value is bounded 
and its importance therefore easily understood. 



Figure 4. The data shown here were computed using lognormal size distributions (w from 1.5 
to 2.5 and r from 0.1 to 87. 1 pm) and Mie theory with ice refractive indices, a) The ratio of 
absorption optical depth to extinction optical depth, R ae (X), versus effective radius for channels 
4 and 5. b) The relationship between R ac (10.8) and Rae(12), which was fit according to: 
Rae(12) = 0.396 -0.287 R ae (10.8) +0.878 R ae ( 10.8) 2 

The Mie-Iognormal model was challenged to reproduce relationships between T 4 and BTD 
measured by AVHRR. For these simulations, it was necessary to vary optical depch to reproduce 
a measured range of brightness temperature. This can be accomplished by fixing number 
concentration and cloud thickness (NAz) and varying r e or by fixing r c and varying NAz (equations 
3, 4, and 10). While results from both approaches can be similar, fixing NAz and varying r e usually 
gives the most convincing reproduction of the measured T 4 versus BTD. This is consistent with the 
understanding of crystal population growth, where particle concentration becomes stable shortly 
after cloud nucleation and particle size continues to increase as a cloud matures. While this 
generalization may not apply to a broad region of cloud, it is probably a good assumption for small 
cloud segments. In general the AVHRR measurements can be explained using the Mie— lognormal 
model, and an example comparison is offered in Figure 5 . 
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T 4 (K) 


Figure 5. Channel 4 brightness 
temperature versus BTD from 
AVHRR cloud measurements at 
5:50 UT on October 7, 1993 over 
the Weddell Sea (roughly 70°S, 
335°E). 855 four km resolution 
AVHRR pixels were used (an area 
of roughly 120 x 120 km). Tbit was 
taken from clear air adjacent to the 
cloud and T c was retrieved from 
AVHRR, as described below. The 
AVHRR data are compared to 
model calculations considering two 
values of NAz and a range of size 
distributions (w =2 and r from 0.1 to 
87.1 pm). 


3. Predictions for PSCs at Equilibrium 


It is useful to estimate optical depths and AVHRR signals expected for typical clouds. Given 
a pressure, temperature, and gas mixing ratios of H 2 0 (Q w ) and/or HNO 3 (Qn), the volume 
densities of stratospheric ice, NAT, and LTA can be calculated from theory by assuming that vapor 
in excess of equilibrium condenses onto cloud particles [e.g., Koop etal., 1997]. This treatment 
is valid because the stability of the stratosphere ensures that vapor loss to cloud growth is localized 
and therefore limited, on the time scale of a cooling event. This approach is invalid in che 
troposphere in principal because a cloud’s vapor supply can be continuous. Infrared extinctions are 
strongly correlated to V, which allows the estimation of optical depths and therefore AVHRR 
signals. Figure 6 a shows an example of the correlation between (3 e (X) and V determined from theory. 
These data were fit according to 

logio(P e (^» = A + B logjo(V) + C logio(V ) 2 ( 11 ) 

where A, B, and C are constants. Constants for (11) are given in Table 2 for ice and NAT and X 
= 10.8 and 12.0 pm. In a different approach, the ratio |3 e (X)/V can be determined versus effective 
radius (Figure 6 b). Either approach is valid, and allows conversion from V to (3 e (X). 
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Figure 6. _The data shown here were computed using lognormal size distributions (w from 1.5 
to 2.5 and r from 0.1 to 87.1 pm) and Mie theory with ice refractive indices, a) Volume versus 
extinction for X = 12.0 pm and a fit to these data (equation 11 ). b) The ratio of extinction over 
volume versus effective radius for X = 10.8 and 12.0 pm. 


Table 2 . Constants for Equation 11 .* 

Wavelength (pm), particle type 

A 

B 

C 

10.8, NAT 

-3.48 

1.11 

-0.054 

10 . 8 , ice 

-3.56 

1.00 

-0.035 

12.0, NAT 

-3.40 

1.04 

-0.045 

12 . 0 , ice 

-3.30 

0.96 

-0.035 

*For P C (X) in km -1 and V in pm 3 cm -3 . 


Volumes corresponding to a range of stratospheric temperatures and pressures were 
computed for upper limit mixing ratios of stratospheric H 2 O (5 ppmv) [Abbas et al., 1996] and 
HNO 3 (15 ppbv) [Roche et al., 1993]. Peak NAT volumes (and therefore optical depths) are roughly 
three orders of magnitude less than peak ice volumes (Figure 7). This difference is due to the fixed 
stoichiometry of NAT so that HNO 3 becomes the limiting reagent for particle growth. Extinctions 
were determined from V, and optical depths were determined for fixed cloud thickness (Figure 7b). 
The optical depths were used with (10) to calculate the change in background sky brightness 
temperature (Lbk(X) = L*(Tbk , X)) when clouds are introduced, Tbk-T m (X), and associated cloudy 
air BTDs (BTD = T 4 - T 5 ) for assumed values of Tbk and Rae(X). The thickest stratospheric NAT 
cloud (a e (X) less than about 10 -2 ) changes the background brightness temperature by less than 0.1 
K and NAT PSCs can therefore be considered invisible to AVHRR. This conclusion holds for a 
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Figure 7. a) Equilibrium volumes for ice and 
NAT determined using 5 ppmv H2O and 15 
ppbv HNO3 (the V. — 1000 pm 2 cnr^ contour 
corresponds to IWP = 0.917 g m -2 , for Az = 
1 km), b) Optical depths determined from 
the volumes in Figure 7a for Az = 1 km, 
using the value of $z(\2)IV from Figure 6b 
for r c = 1 pm for ice. A similar curve was 
used for NAT. The NAT contours were 
terminated 3 K below the ice point 
temperature in Figures 7a and 7b. c) BTDs 
and the difference T bk - T m (12) determined 
for the ice volumes in Figure 7a using (10) 
with Az = 1 km, T bk = 250 K, R ae (X) = 0.6, 
and values of p e (X)/V from Figure 6b for r c = 
1 pm. 
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broad range of conditions (Rae(X), T bk , and T c ), and the same conclusion is drawn for PSCs 
composed of LTA. Considering that the stratospheric conditions in Figure 7 are generous in terms 
of pressure and mixing ratio (in reality HNO3 disappears near the tropopause), this conclusion is 
well founded and the remainder of this study proceeds without consideration of NAT or LTA. The 
thickest ice PSCs change T bk by up to 16 K and produce BTDs of up to 14 K, for the scenario 
presented in Figure 7c. T bk — T m is less sensitive to the shape of the size distribution than BTD, 
and the estimates of BTD in Figure 7c can change notably for using different size distributions, due 
to the dependence of R a c(X) on particle size (Figure 4a). Calculations considering tropospheric 
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conditions produce a broad range of optical depths and BTDs, which include stratospheric values. 
Thus, a threshold optical depth could be determined which precludes a cloud from being 

stratospheric. However, a threshold which precludes a cloud from being tropospheric is more 
elusive. 

It is useful to consider the relationship between cloud optical depth, cloud temperature, and 
brightness temperature. Figure 8 shows example calculations for a single cloud layer in an otherwise 
clear atmosphere. In order for a cloud to radiate at its physical temperature (T m (X) = T c ), its optical 
depth must be greater than about one. Optical depths which produce T m (X) = T c decrease with 

decreasing Tbk (or increasing T c ), although they are often near the upper limit for ice PSC optica] 
depth (Figure 7b). 




Figure 8. Channel 5 cloudy sky brightness temperatures computed from (10) for Rae(12) = 0.6 
and two values of Tbk, contoured for a range of cloud temperatures and optical depths. The 
lower threshold in optical depth where cloudy sky brightness temperature equals physical cloud 

temperature (T m ( 12) = T c ) is indicated. 


4. Cloud Property Retrievals 


4.1. Cloud Temperature and Optical Depth 

Given cloudy and background sky measurements from both channels, it is possible to 
retrieve T c , a c (10.8), and a e ( 12.0) simultaneously. This is accomplished by solving (10) for the 
blackbody cloud radiance L*(T C , X), 
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( 12 ) 


L* (T >0 = ~ I'mW °p( — a »W] 

° 1 - ap[-<r,(AV?„(A)] 

and finding T c through the Planck function. Cloud temperatures can be determined from (12) for 

a range of synthetic optical depths and fixed R ae (X). This step is performed for each channel and 


a solution is taken where both channels yield the same cloud temperature. Solutions are not unique, 
yet the range of possible solutions usually contains the exact solution. It is difficult, however, to 
discern the correct solution from this family of possibilities. One approach is to find the average 


in 


temperature, and use the corresponding optical depths. The retrieval approach is demonstrated i 
Figure 9a, where cloudy sky radiances were simulated for known conditions and (12) was solved 
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Figure 9. a) Optical depths and cloud temperatures which satisfy (12) for channels 4 and 5, 
for cloudy sky radiances simulated using one optical depth pair. The range of possible solutions 
is indicated. One solution is found where a constant temperature intersects the curves for each 
channel. The average solution is taken at the midpoint of the range in T c . b) As in Figure 9a 
except for cloudy sky radiances simulated using five optical depth pairs, and the curves for 
channel 5 are not shown. The exact solutions in Figures 9a and 9b were based on T g = 195 K, 
T bk = 250 K, a s (12.0)/o e (10.8) = 1.5, R ae (12} = 0.75, and R ac (10.8) = 0.82. 


for both channels to determine T c , a e (10.8), and a e (12.0). While an average solution can be a good 
approximation in some cases, it can be biased by over 10 K in others. The assumption is now made 
that horizontal variations in cloud radiance over small distances are due to changes in optical depth, 
while temperature remains fairly constant. This is an important characteristic because for cloud 


radiances corresponding to a range of o e at constant T e the families of solutions tend to overlap 
around the exact value of T c . This behavior is demonstrated in Figure 9b, where the potential 
improvements for considering a range of optical depth are apparent The errors for retrievals based 
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on a single measurement are greatest at low to intermediate optical depths, while T c can be 
accurately determined for optically thick clouds (Figure 9 b). Uncertainties in retrieved T c and a e 
were determined for a range of 7 C and optical depth ratio, R S4 = a c (12.0) / a e (10.8). In general, 
R 54 ranges from two for r e = 0.1 pm to unity for r e = 100 pm. The fundamental uncertainties in 
retrieved T c and o e are systematic (Figure 10) for a chosen scenario of cloud properties. Of course, 
since the cloud properties are unknown for a real measurement, the uncertainties will take on a more 
random nature. 
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Figure 10, Cloud temperatures and extinction optical depths were retrieved from simulated 
radiances based on four cloud optical depths (a e (10.8) = 0.5, 1.0, 1.5, 2.0 and o e (12.0) = 
^ 54 o c( 10 . 8 )) at constant T c . Differences between the retrieved and exact values are shown for 
T c and a e (10.8), for a range of T c and optical depth ratio. The differences in o e (10.8) are an 
average for the four optical depths. The radiance simulations used T b k = 250 K, R ae (12) = 0.75, 
and R ac (10.8) = 0.82. 


4.2. Effective Radius and Ice Water Path 

As illustrated in Figure 6 , the ratio f3 c (X)/V can be a unique function in effective radius. 
Furthermore, it can be shown that the two functions, f(r e , 10.8) = V/p e ( 10 . 8 ) and f(r e , 12.0) = 
WP C ( 12 . 0 ), are monotonic and non-intersecting for roughly 2 < r e < 100 pm. Curves for f(r e , X) 
can be generated using a broad range of size distributions and thus r e . Scaling f(r e , X) by a measured 
optical depth yields a relationship between and r c and VAz, f s (r c , X) = f(r c , X)o e (X) = VAz. Because 
the atmospheric values of VAz and r c will be unique for a given measurement, a broad range of size 
distributions can be searched for one that yields the same r e and VAz from both channels, f s (r c , 10.8) 
= f s (r e , 12 . 0 ), 
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( 13 ) 


1^£ a (12 0) — a no o\ 

^(12.0 ,r e )Az ae(12 - 0) <7e(10 - 8 > 

where the optical depths are from AVHRR. For illustration, f(r e , 10.8) and f(r e , 12.0) were 


calculated using Mie theory with lognormal size distributions. These curves were scaled by realistic 
channel 4 and 5 optical depths, and solutions for r e and VAz are found at their intersection (Figure 
11). The ice water path is easily calculated, IWP = VAz Qj. While size distributions that satisfy (13) 
can vary, the integral quantities tend to fall in a narrower range. Never the less, the curves describing 

fs( r e» ^) tend to have steep intersections, so that solutions become highly sensitive to errors in optical 
depth. 



Figure 11. The ratio of volume over 
extinction computed for lognormal size 
distributions (w=2 and r = 0. 1 to 87. 1 
pm) using Mie theory with ice 
refractive indices. These curves were 
scaled by realistic cloud optical depths, 
Oc(10.8) = 0.2 and ct c ( 12.0) = 0.3, to 
give VAz. The intersection of the scaled 
curves (dotted lines) yields solutions for 
effective radius (8 pm) and VAz (1283 
pm 3 cm -3 km). 


To arrive at the curves in Figure 1 1, the size distribution width was fixed and the median 


radius was varied. This scenario yields one solution. However, a variety of realistic scenarios can 
yield solutions. Sensitivity of r e and VAz to the choice of size distribution scenario was evaluated 
by generating optical depths and corresponding valnes of r c and VAz for a broad range of size 
distributions. These synthetic optical depths were then used with the curves in Figure 1 1 to retrieve 
r e and VAz, for comparison to the known values. For each value of r e> the uncertainty analysis used 
a Monte Carlo approach considering randomly chosen distributions which satisfy (5). Random and 
systematic variability about a known solution were then calculated from this sample set, to 
determine the effects of searching a limited set of size distributions for a solution. The effects of 


10% random errors in the retrieved optical depths were then added to the size distribution errors for 
these calculations. Random uncertainties in r c due to uncer tain ties in the size distribution are less 
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than 30% for the range of r e considered (Figure 12a). When errors in retrieved optical depth are 
added, random uncertainties in r e increase to over 50%, and significant retrieval bias appears for 
r e < 3 pm and r e > 30 pm (Figure 12a). Uncertainties in r e and VAz due to size distribution and optical 
depth errors are summarized in Figure 12b, showing nearly 100% uncertainty in both for some size 
regimes. While r e and VAz cannot be accurately determined, limits can be placed on their possible 
range, as concluded by Parol et al. [1991]. These large uncertainties may be tolerable for some 
applications, and averaging retrieved values over many pixels may provide some reduction in 
uncertainty. 



effective rodius (jim ) 


100 


Figure 12. a) Random and systematic 
uncertainties in r e determined using Monte Carlo 
calculations considering errors in the size 
distribution, n(r) (w varied randomly between 1 
and 3 with the corresponding r determined by 
(5)), and errors in both n(r) and 10% errors in the 
optical depths, b) The nominal values of r c and 
VAz (heavy solid line) and random uncertainties 
(due to both n(r) and optical depth) indicated as 
vertical bars and lighter shading for VAz and 
horizontal bars and darker shading for r c . 


5. Results Using AVHRR 


Cloud properties retrieved from AVHRR were compared to values derived using 
measurements from radiosondes and from POAM II. For this purpose, raw AVHRR data were 
calibrated and geo-registered at 4km horizontal resolution [e.g., Lubin and Harper, 1996]. 

As previously discussed, equilibrium theory allows stratospheric ice clouds to be 
characterized from temperature, pressure, and water vapor mixin g ratio (Q w ). While radiosondes 
provide reliable temperature and pressure measurements, they cannot measure Q w in the 
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stratosphere. Stratospheric water vapor, however, is fairly uniform and has a well known upper 
limit near -5 ppmv [e.g., Abbas et al ., 1996]. Assuming constant Q w , stratospheric radiosonde 
profiles were used to calculate V and (3 e (X) (as in Figure 7), which were integrated vertically to yield 
IWP and o e (X). T c was then computed from a volume weighted mean for all cloud layers. The 
abundance of radiosonde measurements and the straightforward analysis provide a useful approach 
for validating retrievals of T c , o e (X), and IWP from AVHRR ice PSC measurements. An example 
radiosonde analysis is shown in Figure 13, where two values of constant Q w were considered. In 
this example, changing Q w from 5 to 3 ppmv reduces o c (X) and IWP by a factor of four, while T c 
is relatively unaffected. Despite sensitivity to Q w , radiosonde IWP and o e (X) provide reliable upper 
limits, and therefore useful validation data. 
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Figure 13. Radiosonde 
temperatures measured over 
Scoresbysund, Greenland 
(70.5°N, 338.0°E) at 12 UT on 
February 19, 1996. 

Stratospheric frost point 
temperatures (Tj ce ) and cloud 
volumes were computed for 
constant water vapor mixing 
ratios of 3 and 5 ppmv. 
Extinctions were determined 
from V using the curves in 
Figure 6b for r e = 4 pm. T c , 
o c (X), and IWT are listed 
(inset) for the range of Q w 
considered. 


temperature (K) 
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AVHRR cloud temperatures were also compared to solar occultation measurements from 
POAM II [Randall et al . , 1996]. Ice PS Cs are identified in POAM II aerosol extinction profiles 
when the signal saturates in the stratosphere, temperatures are below T; ce , and the termination 
altitude is above the coldest layer in the profile [ Fromm et al., 1999]. When these criteria are met, 
PSC tops are estimated to occur from 1 to 2 km below the profile termination altitude, and PSC 
temperature can be estimated from the corresponding temperature profile. 
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A case study of clouds over Scoresbysund, Greenland beginning on February 19, 1996 was 
used to test the AVHRR retrievals by comparing them with radiosondes and POAM II. AVHRR 
imagery on this day show clouds with brightness temperatures (T 5 ) between 189 and 200 K over 


much of Greenland. A cloudy area within roughly 200 km of Scoresbysund was selected from each 
AVHRR image for the cloud property retrievals. Each AVHRR cloud area was comprised of about 
100 to 300 pixels (an area of roughly 1600 to 4900 km 2 ), and the cloud temperature retrievals used 
all pixels in an area to ensure variability in opdcal depth for reduced retrieval errors (e.g. Figure 9b). 
Values of r c , o e (X), and IWP are presented as averages of all retrievals in a cloud area. Tim* series 


from AVHRR and the radiosondes are compared in Figure 14. These time series depict the passage 
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Figure 14. Time series over 
Scoresbysund, Greenland 
beginning at 00 UT on 
February 19, _1996. AVHRR 
retrievals of T c , a e (X), and 
IWP are compared to values 
derived from radiosonde 
measurements in the 
stratosphere (using Q w = 5 
ppmv; e.g. Figure 13). Error 
bars on the radiosonde values 
correspond to uniform ±2 K 
temperature changes to the 
radiosonde profile. A PSC 
temperature from POAM II 
measurements near 

Scoresbysund is also shown. 
The upper limit for ice PSC 
temperature is indicated in the 
bottom panel (horizontal line, 
193 K). 


of various PSCs and cirrus over Scoresbysund, rather than the evolution of a single starionary cloud. 
The AVHRR cloud temperature at 3 UT on the 19th is too warm to be stratospheric, even when 
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considering a possible 10 K overestimation at low optical depth (e.g.. Figure 9). The other AVHRR 
cloud temperatures are consistent with PSCs, and are in good agreement with coincident POAM 
II and radiosonde data. The AVHRR retrievals of r c , o e (X), and IWP are consistent with PSC 
conditions, yet they do not exclude the possibility of cirrus. Although AVHRR o e (X) and IWP 
retrievals tend to be larger than the radiosonde values, the two data sets agree within the mutual error 
limits. 

Lubin and Harper [1996] suggested that BTDs greater than about 5 K may be characteristic 
of ice PSCs. Their idea was tested in this work using the relationship between measured BTD and 
retrieved cloud temperature from the AVHRR measurements over Scoresbysund. For these 
measurements, stratospheric frost point temperatures (Ti ce ) were always less than about 193 K (e.g. 
Figure 13). Thus the Scoresbysund AVHRR data suggest that ice PSCs are associated with BTDs 
greater than about 5 K (Figure 15), consistent with the suggestion of Lubin and Harper. The use 
of BTD as a PSC identifier may involve more than simply defining a threshold value, and this idea 
is currently the subject of rigorous evaluations based on model calculations and observations. 


Figure 15. Mean BTD versus retrieved 
cloud temperature for AVHRR 
measurements over Scoresbysund on 
February 19, 1996. Mean BTDs were 
computed within each cloud area and 
vertical lines show the range of BTD 
within an area. 


6. Summary 



AVHRR channel 4 and 5 radiances from clear and cloudy skies were evaluated from a 
theoretical basis and techniques were presented for inferring cloud optical depth at X= 10.8 and 12 
pm, average temperature, average effective radius, and ice water path, for ice PSCs or cirrus. The 
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cloud optical depth and temperature retrievals assume that upwelling radiation beneath a cloud layer 
can be characterized from adjacent measurements, and that radiance variations across small cloud 
segments are due to variations in optical depth at constant average cloud temperature. The retrieval 
of ice water path and average effective radius requires signal simulations, which were based on Mie 
theory and lognormal size distributions. While the retrieval assumptions and signal simulations can 
have obvious limitations, they are justifiable for a useful portion of the AVHRR data set. Initial 
assessments show that AVHRR cloudy sky measurements can be explained in terms of the 
Mie-lognormal model. Additionally, cloud temperatures, ice water paths, and optical depths 
retrieved from AVHRR are in good agreement with values derived from coincident radiosonde 
measurements. 

While large values of optical depth, effective radius, and ice water path can preclude the 
possibility of PSCs, small values do not bar the possibility of cirrus. Thus retrieved cloud 
temperatures should provide the most useful information for separating PSCs and cirrus in AVHRR 
data. The retrieval methods described in this paper require a considerable amount of image 
interpretation to determine an accurate underlying background temperature as input to the 
calculations. While this approach is useful for specific case studies, in its current form it would be 
difficult to employ on the entire AVHRR data set. There is compelling evidence, however, that 
could be a straightforward PSC identifier which can be easily applied to the extensive AVHRR 
data set, allowing the development of first-ever PSC climatologies in both hemispheres. 

X 
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